Optical excitations in a one-dimensional Mott insulator 
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The density-matrix renormalization-group (DMRG) method is used to investigate optical exci- 
tations in the Mott insulating phase of a one-dimensional extended Hubbard model. The linear 
optical conductivity is calculated using the dynamical DMRG method and the nature of the low- 
est optically excited states is investigated using a symmetrized DMRG approach. The numerical 
calculations agree perfectly with field-theoretical predictions for a small Mott gap and analytical 
results for a large Mott gap obtained with a strong-coupling analysis. Is is shown that four types 
of optical excitations exist in this Mott insulator: pairs of unbound charge excitations, excitons, 
excitonic strings, and charge-density-wave (CDW) droplets. Each type of excitations dominates the 
low-energy optical spectrum in some region of the interaction parameter space and corresponds to 
distinct spectral features: a continuum starting at the Mott gap (unbound charge excitations), a sin- 
gle peak or several isolated peaks below the Mott gap (excitons and excitonic strings, respectively), 
and a continuum below the Mott gap (CDW droplets). 
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I. INTRODUCTION 

In recent years, various quasi-onCpdjmcnsional mate- 
rials, such as,|eonjugated golymers,BB organic charge- 
transfer salts, &□ Cu oxides, Q and Ni halides,cl have been 
extensively studied because of their unusual optical prop- 
erties and their potential application in modern optical 
technology. In first approximation these materials can 
be described as one-dimensional strongly correlated elec- 
tron systems with half-filled bands. The electron-electron 
interaction dmjes such a system into a Mott insulating 
ground stateoH and dominates low-energy excitations. 
Therefore, the optical properties of one-dimensional Mott 
insulators are currently a topic of great interest. 

Despite many theoretical studies, our knowledge of 
these systems is still fragmentary because of the diffi- 
culties associated with the investigation of strongly cor- 
related systems. For many years, numerical exact di- 
agonalization of small systems was the only method pro- 
viding reliable information on excited states in correlated 
electron systems .Enid Recently, however, the linear optical 
conductivity and exciton properties of one-dimensional 
Mott insulators have been calculated analytically-ji^an. 
infinite system in the limit of a sjxiajl Mott gapE2lll3't3 
and of a large Mott gap.H5E3li3'HE!l Moreover, re- 
cent developments of the density-matrix renormalization 
group (DMRG) methooOEl allows one to calculate ex- 
cited states and dynamical response functions numeri- 
cally in large systems and with an accuracy comparable 
to exact diagonalizationsliiEd 

A paradigm of a one-dimensional Mott insulator is the 
extended Hubbard model with hopping integral f, on-site 
repulsion U, and nearest-neighbor repulsion V at half fill- 
ing. Although this model has been widely studied, its 
properties are still poorly understood in the thermody- 
namic limit. It is knowrcj that the system is a Mott 
insulator in a large region of the parameter space (U, V) 
which is physically relevant for the Coulomb repulsion 
between electrons (U > V > 0). The precise ground- 
state phase diagram has only recently be determined us- 



ing DMRG II In the special case of the Hubbard modelEfl 
(V = 0,U > 0) the ground state is known exactly to be 
a Mott insulatorc3 and the optical conductivity has been 
calculated. til In principle, the optical properties of the 
Mott insulating phase are also knowHJpcJZ.-is_0.at weak 
couplingBB and at strong couplingJiiMBBl-Zl but the 
range of validity of these results is not clear a priori. 
Results for intermediate coupling and close to the Mott 
phase boundary U ~ 2V are scarce and the absence of 
finit e - i j iz Q- Q ffact analysis often hinders their interpreta- 



tion! 

In this paper I present an accurate and comprehen- 
sive investigation of the linear optical conductivity <Jx(uj) 
and the excited states contributing to cti(cj) in the Mott 
insulating phase of the one-dimensional extended Hub- 
bard model at half filling. An efficient symmetrized 
DMRGB *n4_fhe recently developed dynamical DMRG 
(DDMRG)tlfEiJ are used to calculate the optically excited 
states and the linear optical conductivity a\{oj) on large 
lattices. The comparisrin. of numerical resultsJipr ai(u>) 
with field-theoreticaHE^I and strong-couplingiJ predic- 
tions confirms both the great accuracy of DDMRG and 
the wide validity range of both analytical methods. I 
have found that four types of excitations determine the 
optical properties in the Mott insulating phase: pairs of 
unbound charge excitations, excitons, excitonic strings, 
and charge-density-wave (CDW) droplets. I will show 
that each type of excitations dominate the low-energy 
spectrum in a particular region of the parameter space 
(U,V) and exhibits a distinct optical spectrum <7i(w). 

The model, the linear optical conductivity <j\(oj), and 
the relevant symmetries are introduced in detail in the 
next section. In Sec. Ill the numerical methods are 



briefly presented, then the estimation of DMRG trun- 
cation errors and finite-size effects are discussed. In 



Sec. IV I describe the four different types of excitations 
which contribute to the linear optical conductivity and 
the corresponding optical spectra in the various interac- 
tion regimes from the limit of a large Mott gap to the 
limit of a small Mott gap. The final section contains the 
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conclusion. 



II. MODEL 

The one-dimensional extended Hubbard model is de- 
fined by the Hamiltonian 



<Y_ Q+i,<r + c l+1 Ko- 



fi = -*E 

+y^(n / -l)(n ;+1 -l). 



(1) 



It describes electrons with spin a =T, j which can hop 
between neighboring sites. Here Cj , c; !(T are creation and 
annihilation operators for electrons with spin tr at site I, 
ni,a = Ci ^ci^ are the corresponding density operators, 
and hi — hi^ + The hopping integral t > gives 
rise to a a single-electron band of width At. The Coulomb 
repulsion is mimicked by a local Hubbard interaction U , 
and a nearest-neighbor interaction V. The physically 
relevant parameter regime is U > V > 0. The number of 
electrons equals the number of lattice sites N (half-filled 
band). This system is in a Mott insulating phase for 
V < V C (U) « U/2 (Ref.||). Precise values of the Mott 
phase boundary V C (U) are given in Ref. 

Note that the chemical potential is chosen in such 
a way that the Hamiltonian (jl]) explicitly exhibits a 
particle-hole symmetry. This Hamiltonian has two other 
discrete symmetries which are useful for optical excita- 
tion calculations: a spin-flip symmetry and a spatial re- 
flection symmetry (through the lattice center). There- 
fore, each eigenstate has a well-defined parity under 
charge conjugation (P c = ±1) and spin flip (P s = ±1), 
and belongs to one of the two irreducible representations, 
A g or B u , of a one-dimensional lattice reflection symme- 
try group. 

Spectroscopy with electromagnetic radiation is a-.com- 
mon experimental probe of solid-state materials. lj The 
linear (one-photon) optical absorption is proportional to 
the real part o\ of the optical conductivity. For w^O, 
<Ti(u>) is related to the imaginary part of the current- 
current correlation function by 



Im{xM} 



(2) 



For lj > the current-current correlation function is 
given by 



X(« > 0) 
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(3) 



where a is the lattice spacing. Here, |0) is the ground 
state of the Hamiltonian H, \n) are excited states of H, 
and Eq, E n are their respective eigenenergies. Although 
•q = + is infinitesimal, a finite value may be used to 
broaden the resonances at hto — E n — Eq and to reduce 
finite-size effects. J is the current operator 



» iaet v-^ /-+ . 

3 = — E 



(4) 



where — e is the charge of an electron. I set a = e = % = 1 
throughout, and t = 1 is used in figures showing the 
optical conductivity, which means that ai(u>) is given in 
units of e 2 a/h and u in units of t/Ti. 

We note that the current operator is invariant un- 
der the spin-flip transformation but antisymmetric un- 
der charge conjugation and spatial reflection. Therefore, 
if the ground state |0) belongs to the symmetry subspace 
v4+ = (A g ,P c ,P s ), only excited states \n) belonging to 
the symmetry subspace B~ = (B u , — P c , P s ) contribute 
to the optical conductivity. According to selection rules, 
the matrix element (0| J\n) vanishes if \n) belongs to an- 
other symmetry subspace. In this paper, the excitation 
energy E n — Eq of the lowest eigenstate with a non-zero 
matrix element (0| J\n) is called the optical gap E opt . 



III. NUMERICAL METHODS 



DMRGEj'Ej is known to be a very accurate method 
for one-dimensional quantum systems with short-range 
interactions such as the extended Hubbard Hamilto- 
nian (|l]). In this work I use three different DMRG tech- 
niques to calculate ground states, excited states, and dy- 
namic response functions. All three techniques are based 
on the finite system DMRG algorithm. 

First, the usual ground-state DMRG method is used 
to calculate the ground state |0) for a fixed number N a 
of electrons of each spin a. This method provides the 
ground-state energy Eq (iVf , N± ) and allows us to calcu- 
late ground-state expectation values (0|O|0) for various 
operators O, such as static correlation functions. The 
(Mott) gap Em in the single-particle density of states 
of a Mott insulator can also be obtained using this ap- 
proach. At half filling (N a = N/2) the Mott gap is simply 
given by 



E M = 2 [ E q (Nt + 1, iVJ - £fc(iV T , JV X ) 



(5) 



because of the charge-conjugation-aymmetry. 

Second, a symmetrized DMRGlJ technique is used to 
calculate the lowest eigenstates |n) in the B~ symmetry 
sector. This method yields not only the eigenenergies 
E n of the lowest optically excited states (in particular, 
the optical gap E op t), but also allow us to compute ex- 
pectation values (n\0\n) and thus to analyze the nature 
of these states. To optimize the DMRG program per- 
formance, my implementation of the charge-conjugation 
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and spin-flip symmetries differs from the original idea 
presented in Ref. |2^. This is explained in detail in ap- 
pendix. 

Finally, the dynamical DMRG methodM'EJ is used to 
compute the optical conductivity (Q) convolved with a 
Lorentzian distribution of width ?y > 0. Comparisons 
with exact results have shown that DDMRG is a very 
reliable numerical method, which yields spectra with an 
accuracy comparable— to. exact diagonalizations but for 
much larger systemsOEj-Ell 

All DMRG methods have a truncation error which is 
reduced by increasing the number m of density-matrix 
eigenstates kept (for more details, see Refs. [13 and 20). 



Varying m allows one to compute physical quantities (in- 
cluding spectra) for different truncation errors and thus 
to obtain error estimates on these quantities. I have sys- 
tematically used this procedure to estimate the precision 
of my numerical calculations and adjusted the maximal 
number m of density-matrix states to reach a desired ac- 
curacy. The largest number of density-matrix eigenstates 
used in this work is m = 1000. For all numerical results 
presented in this paper DMRG truncation errors are neg- 
ligible. 

All numerical calculations have been performed on lat- 
tices with an even number N of sites using open boundary 
conditions. As we are interested in the properties of the 
Hamiltonian ([!]) in the thermodynamic limit, numerical 
calculations have always been carried out for several sys- 
tem sizes N in order to investigate finite-size effects. The 
largest system size used here is N = 512. If necessary, the 
results have been extrapolated to the infinite system limit 
N — ► oo. To evaluate finite-size effects in a continuous 
spectrum one has tOpCompute it for different sizes while 
keeping -qN = const.Eil In this work 7]N — 12.8t is used. 
For all numerical results presented in this paper finite-size 
effects (including chain-end effects) are negligible unless 
discussed explicitly. For spectra this means that finite- 
size effects are completely hidden by the broadening rj. 
More precisely, DDMRG results for finite N = 12.8t/rj 
are not distinguishable from the corresponding infinite- 
system spectra convolved with a Lorentzian distribution 
of width r\ (see the discussion in Ref. ^l]) . 



IV. RESULTS 

To facilitate the comparison with analytical results, 
the discussion of optical excitations in the Mott insulat- 
ing phase is divided in three subsections: the limit of a 
large Mott gap, the regime of finite Mott gaps, and the 
limit of a small Mott gap. Note, however, that the Mott 
gap Em just fixes the energy scale; the minimal energy 
required to create a charge excitation is Em/ '2 but op- 
tical excitations do not differ qualitatively as Em varies 
if everything else is kept constant. In all cases, there 
is spin-charge separation and the spin sector is gapless. 
Elementary excitations in the charge sector are spinless 
bosons in the lower and upper Hubbard bands. Optical 



excitations are always made of an even number (> 2) 
of elementary excitations with opposite charges (to pre- 
serve charge neutrality). The different types of optical 
excitations and optical spectra found in the model (Q) 
result from the residual interactions (essentially the non- 
local part of the Coulomb repulsion, here V) between the 
elementary charge excitations. 



A. Limit of a large Mott gap 

In the strong-coupling limit U 3> t, the properties of 
the model |l]) in the Mott insulating phase can be de- 
scribed using simple concepts. In the ground state double 
occupation is prohibited and there is exactly one electron 
on each site. Elementary charge excitations can be rep- 
resented as an empty site (holon in the lower Hubbard 
band) or a doubly occupied site (doublon in the upper 
Hubbard band). The minimal energy required to create 
a holon or a doublon is Em/2 = U/2 — 0(t) ^> t. Optical 
excitations always consist of an equal number of holons 
and doublons to conserve the total charge. The ionicity 
of excited states is defined as the change in the number 
of doubly occupied sites with respect to the ground state 



(n\N d \n) - (0\N d \0) , 



(6) 



where Nd — Yli ^i.t^M an< ^ denote the ground 

state and excited states, respectively. Thus, I n is a mea- 
sure of the number of doublons (or equivalently of holons) 
created by an excitation. Depending on the strength of 
the nearest- neighbor interaction parameter V, the low- 
energy optical excitations are made of a single doublon- 
holon pair (I n — 1) or are collective excitations of several 
such pairs (7„ > 1). Note that /„ is also equal to the 
derivative of the excitation energy E n — Eg with respect 
to U. 



Single holon- doublon pair 

For V < U/3 + 0(t), optical excitations consist of a sin- 
gle holon-doublon pair and-the optical properties, which 
can be calculated exactlyJiJo depend only on the param- 
eters V and t. For < V < 2t, holon and doublon are 
independent. A schematic representation of this state is 
shown in Fig. [I]. This pair of free charge excitations gives 
rise to a continuous band in the optical spectrum ai{u>). 
The band starts at the Mott gap Em — U ~ At and has 
a width of 8t. As there is no optical excitation with a 
lower energy than Em, the Mott gap is also the optical 
gap E opt . The optical spectra for V = and V = 2t 
are shown in Fig. |2| with a broadening rj/t = 0.1. At the 
conductivity threshold cxi(w) vanishes as \Jio — Em for 

V < 2t but diverges as — Em in the special case 

V = 2t. The optical conductivity also has a small peak at 
lu = U — V with 1 % of the spectral weight Jl3 This peak 
is visible inside the band for V = and V = 2t in Fig. |^. 
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FIG. 1: Schematic representation of (a) an unbound holon- 
doublon pair, (b) an exciton, and (c) a biexciton in the strong- 
coupling limit U — 2V ^> t. 



It corresponds to .-a, bound state made of dispersionless 
charge excitationaiil and can be seen as a localized exci- 
ton with size £ = 1. 

For V > 2t, there is also a continuous band start- 
ing at lu = Em = U — At due to independent holon- 
doublon pairs, but the lowest optical excitation is now 
an exciton (i.e., a bound holon-doublon pair) with an 
energy w cxc = E opt = U - V - At 2 /V (Ref. [TJ). The 
term — At 2 /V is the kinetic energy lowering due to the 
exciton center-of-mass motion. Therefore, in the strong- 
coupling limit (U 3> t) the exciton binding energy is 
E b = E M - E opt = V-At + At 2 /V {V > 2t). This bind- 
ing energy significantly differs from the .incorrect result 
Eb = V often reported in the literaturejHj which is (ap- 
proximately) valid only under the additional condition 
V t. In the optical spectrum <j\ (u>) the exciton gen- 
erates an isolated (5-peak at o; GXC below the band onset. 
For V 3> t, an exciton is essentially the nearest-neighbor 
holon-doublon pair shown in Fig. [TJ. This is exactly the 
state generated by the current operator (|]) applied to 
the ground state with one electron on each site. Thus 
the spectral weight is concentrated in the excitonic peak 
for F>i. This strong excitonic peak is already clearly 
visible for V = bt in Fig. g. For finite V/t, however, 
there is a finite probability of finding holon and doublon 
at a distance m > 1. This probability can be calculated 
exactlyE3 



C(m) = C(l - S m , )e 



(7) 



with k = 2 ]n(V/2t) and a normalization constant C . The 
exciton size (the average holon-doublon distance) is then 



t(V > 2t) 



V 2 



V 2 - At 2 



(8) 



and decreases as V increases. Correspondingly, one 
observes a progressive transfer of spectral weight from 
the band above Em to the excitonic peak at lo cxc — 
U - V - At 2 /V as V increases. Note that, as for V < 2t, 
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FIG. 2: Reduced optical conductivity wcri(w) in the limit of 
a large Mott gap (U > t) calculated with DDMRG for three 
different values of V using r\jt = 0.1 (N = 128 sites). 



there is a small peak corresponding to a localized exciton 
at uj = U — V in the optical conductivity. This peak lies 
in the band for V < At but is situated between the band 
and the strong excitonic peak for V > At. In Fig. it 
can be seen for V — 5t as a small bump at the foot of 
the strong peak. 

Figure p shows optical spectra calculated with the 
DDMRG method on a 128-site lattice. On the scale 
of this figure, the DDMRG spectra are indistineuisliable 
from the analytical results for an infinite system.EjIiJ (For 
an expanded view showing small deviations, see Fig. 1 of 
Ref. [13].) This perfect agreement confirms the validity of 
the strong-coupling calculations done in Ref. [TJ]. More- 
over, it confirms once more that DDMRGpGan accurately 
reproduce infinite-system optical spectra.Eil 



Collective excitations 

As we have just seen, the creation energy of an excited 
state with ionicity /„ = 1 is either U — 0(t) for an un- 
bound holon-doublon pair or U — V — 0(t 2 /V) for an 
exciton. Once a first excitation has been created, how- 
ever, the creation of a second holon-doublon pair bound 
to the first excitation requires only an energy U — 2V. 
Therefore, when V becomes large enough, the lowest op- 
tical excitations are bound states of n BY i- excitons, called 



excitonic strings S3c3 A biexciton (n cxc = 2) is shown in 
Fig. [I] as an illustration (see also Ref. [)]). 

The excitation energy of a n CX c-exciton string is 

E{n cxc ) = U-V + (n cxc -l)(U-2V)-0(t 2 /V) . (9) 

Here the correction of order t 2 /V corresponds to the ki- 
netic energy lowering due to the center-of-mass motion. 
An n cxc -exciton string (n cxc > 2) appears in the low- 
energy excitation spectrum, around or below the onset 
Em = U — 4i of the band of free holons and doublons, if 
E{n exc ) < E M or 



V> 



1 



U + 0(t) . 



(10) 
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Thus, the biexciton becomes a low-energy excitation for 
V > U/3 and longer excitonic strings (n exc > 3) for 
larger V. The = 1 corresponds to the usual 

exciton, which appears in the low-energy spectrum as 
soon as V > 2t as discussed above. As an n oxc -exciton 
string is made of n cxc doublons and holons bound to- 
gether, it is a neutral excitation, its ionicity is I = n cxc , 
and its length is 2n cxc — 1 in units of the lattice con- 
stant. Excitonic strings have been observed in the non- 
linear optical spectrum of quasi-one-dimensional neutral 
mixed-stack charge-transfer solids and are known to con- 
tribute to the non-linear optical conductivity of models 
such as the extended Hubbard model ([!]) .BO Naively, one 
does not expect excitonic strings with / = n cxc > 2 to 
contribute to the linear optical spectrum ai(u>). In the 
limit U/t 3> t the current operator (^) creates at most 
one holon-doublon pair and thus in Eq. (^) the matrix 
elements (n\ J|0) between an excited state \n) and the 
ground state |0) must vanish if the ionicity (o) is larger 
than 1. Yet, we will see in the next section that excitonic 
strings with nexc > 2 are visible in the linear optical con- 
ductivity of the extended Hubbard model for large but 
finite couplings U/t and V/t. The reason is that for any 
finite t there are quantum charge fluctuations (virtual 
holon-doublon pairs) in all eigenstates of the Hamilto- 
nian (|l|) which leads to small but finite matrix elements 
(n\ J|0) even if the average ionicity /„ of an excitation \n) 
exceeds one, 



(n\J\0) 



t 



U-2V 



(In-l) 



(11) 



Thus, for V > U/3 3> t the low-energy optical spectrum 
consists of a strong excitonic peak at lo — E op t — U — 
V - 0(t 2 /V) followed by several [n exc = 2, 3, . . . < (U - 
V)/(U—2V)] weaker isolated 5-peaks with exponentially 
decreasing spectral weight atws; E opt + (n cxc — l)(U — 
2V) . All these peaks appears below (or about) the onset 
of a weak continuum due to free holons and doublons at 
lo = E M - 

As long as U — 2V ^> t, excitonic strings retain a well- 
defined size represented by an integer number n oxc be- 
cause the kinetic energy lowering ~ t due to size fluctu- 
ations is much smaller than the energy cost ~ (U — 2V) 
per exciton in the string. Close to the phase boundary 
(U « 21/Lbgtween the CDW ground state and the Mott 
insulator,E3E3 however, size fluctuations become advan- 
tageous. Thus, for U — 2V < t low-energy excitations 
(of the Mott insulator) are CDW droplets, which can be 
understood as superpositions of excitonic strings of every 
size, 



= cx\n cxc = 1) + c 2 |n cxc = 2) + c 3 |n oxc = 3) + 

with a broad distribution of coefficients c n (%2 n 
1). For U — 2V — » + the distribution becomes flat, i.e., 
\ c n\ 2 — * const. (For comparison, an n GXC -exciton string 
can be described by the above state with |c„| 2 w 1 for 
n = n cxc and |c„| 2 <C 1 for n ^ n GXC .) CDW droplets in 



the Mott insulating phase are the analogue to the 3PW 
droplets in the CDW insulator discussed by Hirsch.c3 As 
excitonic strings, these CDW droplets are neutral ex- 
citations, but one can generalize the concept to CDW 
droplets carrying charges (see below). The average size 
rcDw °f a CDW droplet is related to its ionicity by 
rcDW = 2/ = 2 J2i *l c i| 2 - Its excitation energy is 

E(r CDW ) = U~V+ 1 ^(U~2V)-st , (13) 

where s > and —st represents the kinetic energy low- 
ering due to droplet size fluctuations and center-of-mass 
motion (s as 4 for U — 2V — ► + ). Contrary to excitonic 
strings, the ionicity / of a CDW droplet is not a inte- 
ger number but can take any value > 1. Therefore, for 
U — 2V <C t there is a band of CDW droplet excitations 
starting at U — V — st. Moreover, the matrix element 
(n| J|0) for a CDW droplet \n) is essentially given by the 
overlap c\ with the single exciton state in Eq. (|l^). Thus 
in this regime one expects that the CDW droplets give 
rise to a band in the optical spectrum o~\{io) starting at 
to = E opt = U —V —st. This band lies below the Mott gap 
Em ■ It should be noted that the Mott gap is determined 
by the excitation energy of unbound holons and doublons, 
E M w U - 4i, as long as U - 2V > t. For U - 2V < t, 
however, CDW droplets carrying a charge ±e [a CDW 
droplet (|l2| ) bound to an extra holon or doublon] have 
a lower energy than a bare holon or doublon and reduce 
the gap for charge excitations (||) to Em — U — s't with 
s' w 8 for V -> U/2. 

In summary, in the limit of a large Mott gap (U 
t) there are four distinct regimes corresponding to four 
types of excitations in the low-energy optical spectrum: 
independent charge excitations for V < 2t, excitons for 
V > 2t but V < U/3, excitonic strings for V > U/3 > 
t but U - 2V > t, and CDW droplets for U - 2V < 
t. We will see that these excitations are also found in 
the optical spectrum of the extended Hubbard model for 
finite interaction strengths and Mott gaps. 



B. Regime of finite Mott gaps 

For finite coupling parameters U and V the low-energy 
optical properties of the extended Hubbard model (|l|) can 
be calculated using the ground-state and symmetrized 
DMRG methods presented in Sec. III. Figure || shows the 
Mott gap Em and the optical gap E opt as a function of V 
for three values of U. Both gaps increase monotonically 
with U but decrease with increasing nearest-neighbor in- 
teraction y .EJEZl For all values of U the optical gap equals 
the Mott gap (in the thermodynamic limit) as long as 

M ■ 



(12) V < 2t but for larger V, E opt becomes smaller than E 

This suggests that for all U > the low-energy excita- 
tions are unbouiidr&r V < 2t and bound for V > 2t as in 
the U 3> t limit Obviously, the condition V > 2t can 
be realized only for relatively strong coupling (U ^ it) 
because the Mott insulating phase exists only for V up 
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FIG. 3: Mott gap Em (upper line) and optical gap i? op t 
(lower line) versus V for U/t = 40 (dot-dashed), 6 (solid), 
and 2 (dashed). For U = 2t, E opt = E M - 



to V c ~ U/2. The Mott gaps in Fig. || are initially almost 
constant as V increases then diminish significantly close 
to the phase boundary V c . This agrees with the strong- 
coupling analysis in the previous section, which suggests 
that Em is essentially independent of V for V c — V ^> t 
but is reduced by quantity oc t as V approaches the crit- 
ical value V c . Note that on the critical line between the 
CDW and Mott insulating phases, both gaps vanish for 
U < 3i while the Mott gap clearly remains finite for 
stronger coupling (U > At). 

To determine the nature of the low-energy optical ex- 
citations I have calculated their ionicity (g). Figure |] 
shows the ionicity I\ of the first optically excited state 
(the 1B~ state) as a function of V for three values of 
U. In the half-filled Hubbard model (V = 0) h in- 
creases monotonically from at U = to 1 for U — > oo, 
reflecting the increasingly ionic nature of the elemen- 
tary charge excitations in the lower and upper Hubbard 
bands. The ionicity increases slowly with V and remains 
below or close to 1 for most couplings (U, V), which con- 
firms that the corresponding optical excitations are made 
of a single pair of elementary charge excitations. In the 
regime U « 2V > At, however, one observes a rapid 
but continuous increase of I\ to values larger than 2 as 

V — > V c ~ U/2. This shows that the lowest optical ex- 
citation has become a CDW droplet. Looking at higher 
optical excitations, one finds the same qualitative behav- 
ior of the ionicity as a function of U and V. Additionally, 
one observes the formation of excitonic string with inte- 
ger I = n cxc > 2. 

To determine whether the pair of elementary charge 
excitations form bound (exciton) or unbound states, one 
can calculate the average distance between . both exci- 
tations using an exciton correlation function.tSHa For 

V < 2t I have found that this average distance always di- 
verges with increasing system size N. This result defini- 
tively confirms that in this regime an optical excitation 
is a pair of independent charge excitations, in agreement 
with the strong-coupling analysis. For V > 2t, the av- 
erage distance tends to a finite value for N — > oo as 
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FIG. 4: Ionicity 7i of the first optically excited state (1B U 
as a function of V for three different values of U. 



expected for an exciton. The exciton size £ determined 
with this procedure is shown in Fig. ^ as a function of 
V for two finite values of U. The exciton size in the 
limit U 3> t, Eq. (jg), is also plotted for comparison. The 
size £ increases and diverges as V tends to 2t, showing 
the unbinding of the exciton at V = 2t. Note that for 
U = 40f [ -t|ie 1 pKes measured with the exciton correlation 
functiont3E3E3 agree perfectly with Eq. (||). 

One can gain some knowledge about the nature of op- 
tically excited states by looking at the scaling of their 
spectral weight 



tt |(n|J|0)p 



E 



Na E n 



(14) 



[see Eqs. (|) and @] with the system size N. For V < 2t 
I have found that the spectral weight of low-lying excita- 
tions vanishes for N — > oo, which is expected for states 
belonging to a continuum. (As there is an infinite num- 
ber of states in a continuous band, the spectral weight of 
each state must go to zero as N — > oo, so that the to- 
tal weight in any finite frequency interval remains finite.) 
For V > 2t (but outside the CDW droplet regime), I 
have found that the optical weight W\ of the 1B~ state 
tends to a finite value in the thermodynamic limit. This 
corresponds to a <5-peak (with total weight > Wi) in the 
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FIG. 5: Exciton size as a function of V for U = 40t (circles) 
and U — 8t (squares). The solid line is the U 3> t result, 
Eq. M. The dashed line is just a guide for the eyes. 
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FIG. 6: Reduced optical conductivity ujai(uj) for U = 40* 
and four different values of V calculated using DDMRG with 
lj/t = 0.4 (N = 32 sites). 

optical spectrum as expected for an exciton or an exci- 
tonic string. In the CDW droplet regime, finite-size ef- 
fects become large and complex and, in most cases, it has 
not been possible to determine the scaling of the spectral 
weights (|l4|). 

Optical spectra 



10' 




FIG. 7: Reduced optical conductivity ujaxiij) for U = 40* 
and (a) V = 16* and (b) V = 19.97* calculated using DDMRG 
with r\jt = 0.4 (A^ = 32 sites). Vertical lines indicate the Mott 
gap E M - 



The above analysis shows that low-energy optical ex- 
citations in the regime of finite Mott gaps are identical 
to those found in the limit of a large Mott gap and can 
be interpreted using the sim ple the ory developed for the 
strong-coupling limit in Sec. IV A . Turning next to the 
optical spectrum I have calculated <j\{uj) for various pa- 
rameters 40* > U > 3f and U/2 > V > using the 
DDMRG method. I have found that the optical spectra 
of systems with a finite Mott gap closely resemble those 
observed in the limit of a large Mott gap. 

As a first example, the optical conductivity <Ji(u>) is 
shown in Fig. || for U — 40* and several values of V 
representing the four different regimes: free charge ex- 
citations (V = 0), excitons (V = 5*), excitonic strings 
(V = 16i), and CDW droplets (V = 19.97i). For 
V = there is a single continuous band starting at 
u) = E opt = Em = 36.14*. At the band edges the op- 
tical conductivity vanishes as \Jio — Em as discussed in 
Ref. 0. 

For V = 5t a strong excitonic <5-peak appears at 
w cxc = E opt = 34.39* below the Mott gap E M = 36.13*. 
The exciton has a size £ w 1.2 in perfect agreement with 
Eq. (^). There is also a weak continuous band of free 
charge excitations above Em, which is only visible as a 
high-frequency tail of the exciton peak in Fig. [j| The 
gap between the excitonic peak and the band is not vis- 
ible in Fig. U because of the large broadening rj/t = 0.4 
used here, but it can be checked with a scaling analy- 
sis for rj — > (N — > oo) as discussed in Ref. ^l| The 
only qualitative difference between the present result for 
U = 40t and the corresponding result in the limit U ^> t 



(see Fig. 1 in Ref. |l3|) is the absence of the weak peak 
associated with a localized exciton at u> — U — V. Nev- 
ertheless, this weak peak is not an artifact of the strong- 
coupling limit because its existence has been confirmed 
in the Hubbard model (V = 0) down to U = it (Ref. 0). 
The finite spectral weight carried by the localized exciton 
originates from a ground-state-dimer-dimer correlation 
of the spin degrees of freedomJiil In the strong-coupling 
limit (U - 2V > t) of the extended Hubbard model (§), 
the effective exchange coupling between nearest-neighbor 
spins depends on the occupation of the neighboring sites 
if V ^ 0. Thus the effective spin Hamiltonian is not the 
one-dimensional Heisenberg model, in general. Only for 
V = or in the limit U ^> V, the effective spin Hamil- 
tonian reduces to the Heisenberg model with a constant 
exchange coupling J = At 2 /U. In this case, the ground 
state has the relevant spin dimer-dimer correlations and 
the localized exciton carried a finite optical weight as ex- 
plained in detail in Ref. 14. For finite U and V, however, 



the spin dimer-dimer correlation is presumably destroyed 
by the fluctuations of the spin exchange coupling and 
thus the optical weight of the localized exciton vanishes. 

For V = 16* the condition V > U/3 is satisfied and 
excitonic strings appear in the optical spectrum below or 
around the Mott gap Em = 35.58*. As seen in Fig. ||, 
most of the spectral weight is concentrated in the exciton 
of size £ = 1.0 at cj cxc = E opt = 23.53* mU-V. The 
biexciton at ui = 32.35* w 2U — 3V is barely visible in 
Fig. |^. The optical conductivity o\ [uS) is again shown 
in Fig. |7](a) on a logarithmic scale. The isolated peaks 
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FIG. 8: Optical conductivity ai(to) for U = 8t and four 
different values of V calculated with DDMRG using rj/t = 0.1 
(N = 128 sites). 



associated with both excitations are now clearly visible. 
The measured ionicity (||) is / = 1.1 and I = 2.2 for the 
exciton and the biexciton, respectively. In Fig. ^(a) the 
remnant of the continuous band of free charge excitations 
and the triexciton (at u> « 3U — 5V = 402) are also visible 
in the interval u> = 36 — 442 above the Mott gap. 

For V — 19.972; « U/2 the optical conductivity spec- 
trum is radically different. The excitonic strings collapse 
into a band of CDW droplets with varying sizes. For in- 
stance, the 1B~ state is a droplet of size j-cdw = 8.9 with 
an energy E opt s» 15.5£. These CDW droplets give rise to 
a broad band in the optical conductivity spectrum shown 
in Fig. |(| The onset of this band is well below the Mott 
gap Em = 30.832. On the logarithmic scale of Fig. 0(b), 
one sees that, in this particular case, the entire optical 
weight seems to be below Em (for 77 — > 0) . The appear- 
ance of a band below the Mott gap is also visible in the 
current-current correlations for U = 12t and V = 62 pre- 
sented in Ref. ^ but the optical spectrum in the regime 
U « 2V is not interpreted correctly in that work. 

As a second example and to illustrate the finite-size- 
scaling analysis I have carried out for dynamical spec- 
tra, I discuss the optical conductivity o\{w) for U = 82. 
Figure || shows the evolution of the optical conductiv- 
ity for increasing nearest-neighbor repulsion V. For 
V = t and V = 2t, the spectrum contains a single con- 
tinuous band due to free charge excitations starting at 
E opt = E M = 4.672 and 4.532, respectively. For V = 3t 
the spectrum consists of a strong peak corresponding to 
an exciton of size £ = 3.2 and energy u> cxc = E opt = 3.862, 
and of a weak band above the Mott gap Em = 4.102. 
This band is visible in Fig. || as the high-frequency tail of 
the excitonic peak. For V = 42, CDW droplets of vary- 
ing sizes dominate the optical spectrum. For instance, 
the 1B~ state is a droplet of size tcdw = 5.6 with an 
excitation energy E opt = 1.552 lower than the Mott gap 
Em = 2.292. There is no intermediate regime with well- 
defined excitonic strings for this value of U. 

The precise shape of 01 (u;) cannot be determined from 
the sole results shown in Fig. because of the finite res- 
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FIG. 9: Scaling analysis of the maximum in DDMRG opti- 
cal spectra cri(w) for U = 8t and three different values of V 
(r/N — 12. Si), (a) Position oj max of the maximum (open sym- 
bols) as a function of the inverse system size. Lines are linear 
fits to these data. Solid symbols show the optical gaps E opt 
calculated with symmetrized DMRG for TV — > 00. For V = 32, 
the Mott gap Em is indicated by a triangle. For V/t — 1 and 
2, Em = E pt. (b) Maximum of ai(u>) as a function of rj/t. 
Lines are fits to the numerical data: the dashed line corre- 
sponds to r? -1 , the dot-dashed to rj~ 1 ^ 2 , and the solid line is 
a linear fit in r\. 



olution and system size used, rj/t = 12.8//V = 0.1. To 
determine the properties of u\{uj) with maximal resolu- 
tion (77 — > 0) in the thermodynamic limit (N — > 00), one 
can perform a scaling analysis with r]N = const, as ex- 
plained in Ref. [H]. (Here I have used nN = 12.82.) The 
scaling analysis of the optical conductivity cti(cj) calcu- 
lated with DDMRG always yields results which are qual- 
itatively and quantitatively consistent with the proper- 
ties of low-lying optical excitations determined using the 
ground-state and symmetrized DMRG methods. For in- 
stance, <j\ (to) vanishes for all u> < E opt and there is a 
continuous band for u> > E opt or a (5-peak at u) = E opt in 
the limit 77 ~ l/N —> 0. 

The scaling analysis of the conductivity maximum 
Cmax = cri(w m ax) in DDMRG spectra is illustrated in 
Fig. ^ for the same interaction parameters as in Fig. ||. 
For V — 2, w max tends to a value (5.152) larger than the 
optical gap E opt = 4.672 calculated with symmetrized 
DMRG for N — > 00 [see Fig. ^|(a)] while er max tends to 
a finite value for 77 — > [see Fig. [|(b)]. Moreover, the 
derivative of <j\ (u>) has a maximum that diverges as l/y/rj 
for r\ — > at lu — E opt . These features correspond to a 
continuum that vanishes as ^Jio — E opt at the conductiv- 
ity threshold and goes through a maximum just above 
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the optical gap at 
V = 2t, uj m ., 
E opt for N 



lAE opt (see Ref. |21|). For 



tends to the same value as the optical gap 
* oo [see Fig. Ej](a)] and er max diverges as 
for rj — > [see Fig. g(b)]. These features corre- 
spond to a continuum thatdiverges as 1/ yfuj — E op t at 
the conductivity threshold. Ell (Note that for all values of 
U investigated <7i(w) displays this divergence at V = 2t.) 
Therefore, the features of the optical spectrum for V < 2t 
are similar to those found in the strong-coupling limit. 
For V = 3i, w max tends for N — > oo to the same value 
as the optical gap E op t (which is smaller than the Mott 
gap in the thermodynamic limit) [see Fig. |^(a)] and <x max 
diverges as 1/r} for rj — > [see Fig. ||(b)]. These features 
correspond to a 5-peak at w = -Eopt- Moreover, oi(w) 
vanishes between E opt and but remains finite above 
Em in the limit r\ ~ 1/iV — > 0. Therefore, in the thermo- 
dynamic limit the spectrum for V = 3t (shown in Fig. |§| 
for N = 128 sites) consists of an excitonic S-peak sepa- 
rated from the band of independent charge excitations as 
in the strong-coupling limit. In the CDW droplet regime 
(i.e., close to the critical line V c » U/2 separating Mott 
and CDW phases) finite-size effects are more complicated 
and larger than in the other regimes. As a consequence, 
for V — At it has not been possible to perform a con- 
clusive analysis with the largest system sizes (N = 256) 
available. It seems that the low-energy spectrum con- 
tains a 5-peak at E opt and a band starting immediately 
above E opt , both due to CDW droplets. 

In summary, I have found that the optical properties 
for all finite Mott gaps (i.e., for all U > 0, V > in the 
Mott phase) are qualitatively similar to th ose calculated 
in the limit of a large Mott gap (Sec. [IV A| , Refs. [b| 

(i) For V < 2t (that is the only possible case for 
U < 4t), independent charge excitations give rise to a 
continuous band starting at the Mott gap Em, which 
is equal to the optical gap. The band width is typically 
~ 8t. For V < 2t, <Ji(ui) vanishes smoothly at the thresh- 
old E opt , typically as */w — E op t ■ At V = 2t, <7i(u>) di- 

-£o P t->0+. 

-0(t) andV" < U/3+0(t) 



E op t for ui 



verges as 1/ */a 

(ii) For V > 2t but U>2V^t 
(this is possible only for U > At), the optical spectrum 
consists of an excitonic 5-peak below the Mott gap and 
a band due to free charge excitations above Em- Most 
of the optical weight is in the excitonic peak for V > it. 

(iii) If U is large enough (U > 12t) excitonic strings 
appear in the low-energy spectrum for V > U/3 + 0(t) 
but U > 2V+0(t). They generate isolated <5-peaks below 
the Mott gap Em in the optical conductivity u\ (u>) with a 
separation between peaks of Aco w U—2V. The first peak 
is an exciton and contains most of the spectral weight. A 
very weak band due to free charge excitations still exists 
above Em- 

(iv) Close to the boundary V c ~ U/2 between the Mott 
and CDW phases, if V exceeds 2t, CDW droplets dom- 
inate the low-energy spectrum and give rise to a broad 
band (including sharp peaks) starting below the gap for 
charge excitations (^). 



TABLE I: Mott gap Em, optical gap E op t, ionicity 7i of the 
first optically excited state 1B~ , and the corresponding field- 
theory interaction parameter (3 2 (see text) for several values 
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In the limit of a small Mott gap (Em <S t) the co- 
herence length ~ At/ Em becomes very large and it is 
not possible to carry out numerical simulations on lat- 
tices large enough (N 3> At/ Em) to determine the opti- 
cal spectrum with confidence. Fortunately, in this limit 
field-theoretical methods provide generic results for the 
low-energ y . puti cal spectrum of a one-dimensional Mott 
insulatorOl!3ll3 Field-theoretical results are applicable 
to lattice models such as Eq. ([!]) for gaps up to Em ^5 t, 
which makes possible a direct quantitative-comparison of 
field theory and DDMRG calculations MBEI 

In the field-theoretical approach, elementary charge ex- 
citations are holons (in the lower Hubbard band) and 
anti-holons (in the upper Hubbard band). Optical exci- 
tations are made of a equal number of holons and anti- 
holons. Assuming that the low-energy excitations consist 
of one holon-antiholon pair, the optical conductivity is 



FT 



M = A Sp{uj/Em) 



(15) 



where Sp(x) is a known function depending on the field- 
theory interaction parameter < /3 < 1, and A is 
a unknown constant which sets the conductivity scale. 
Strictly speaking, this result is exact only for ui < 2Em 
and 1 > 1/3 but it has been found by comparison with 
DMRG results that corrections for oj > 2Em are usually 
negligible For 1 > 1/2, <7f T (w) describes a single 
continuous band starting at E op t = Em, which is due to 
independent holons and anti-holons. The optical conduc- 
tivity vanishes smoothly as J~w — E op t for fi 1 > 1/2 and 
diverges as 1 / y/ui — E op t for 1 — 1 /2 at the conductiv- 
ity threshold. For 1/3 < (3 2 < 1/2, there is a <5-peak 
at lo = E pt < Em in addition of the band starting at 
lu = Em- The <5-peak is due to a bound holon-antiholon 
pair (exciton). For f3 2 < 1/3 additional excitons and 
excitonic strings (made of several holon-antiholon pairs) 
appear in the spectrum. Therefore, field-theoretical pre- 
dictions for the optical conductivity of a one-dimensional 
Mott insulator are qualitatively similar to what we have 
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FIG. 10: Optical conductivity ai(iv) calculated with 
DDMRG using rj/t = 0.1 (N = 128 sites) for U = 3t, V = 
(dashed), U = 3.5*, V = 1.4* (dot-dashed), 17 = 4.15*, V = 2* 
(dotted), and U = 6t,V = 3.155* (solid). 



found in the extended Hubbard model (|l]) using a strong- 
coupling analysis and DDMRG simulations. 

The field-theory parameters Em, 13, and A must be es- 
timated numerically by comparison with DMRG results 
because one does not know their relations to the lattice 
model parameters U, V, t. I have first determined sev- 
eral couplings (U, V) which yield approximately the same 
Mott gap Em /* ~ 0.6 — 0.7. These couplings are listed in 
Table | with the corresponding Mott and optical gaps cal- 
culated using the ground-state and symmetrized DMRG 
methods. Then I have calculated the optical conduc- 
tivity <Ji(u) for these parameters using DDMRG. Some 
results are shown in Fig. |l0|. Note the progressive dis- 
placement of spectral weight to lower energy as V in- 
creases although the Mott gap remains almost constant 
(see Table |). To determine the parameters (3 and A 
one can now compare the field-theoretical spectra with 
the DDMRG data. [<7f T (w) has to be convolved with a 
Lorentzian distribution of appropriate width r\ to make 
a direct comparison. Ell] This procedure yields the param- 
eter (3 2 listed in Table Q. As expected the boundaries 
between independent charge excitations and excitons in 
the field theory and in the lattice model ([!]) coincide: 
1 = 1/2 corresponds to V = 2*. 

For V < 2t (corresponding to (3 2 > 1/2) one can find 
parameters so that trf T (a;) perfectly fits the numerical 
data over a wide frequency range. For instance, in Fig. |ll| 
no difference is visible between the DDMRG spectrum 
for U = 4.15* and V = 2* and the fitted field-theoretical 
spectrum up to ui = 2* ss 3Em- For V > 2*, however, 
discrepancies between DDMRG and field-theory results 
appear and grow progressively stronger as V increases. 
It is no longer possible to find parameters j3 and A to 
reproduce the DDMRG spectra over a significant fre- 
quency range above ui — E opt . Instead (3 and A are set 
by the optical gap and the total spectral weight. This 
yields the values of (3 shown in parenthesis in Table |[ As 
an example, one see in Fig. [ll] that the field-theoretical 
spectrum differs significantly from the DDMRG result 
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FIG. 11: Comparison of optical spectra calculated with 
DDMRG (N = 128 sites) and field theory (N = oo) for 
r]/t = 0.1. The dashed line is the DDMRG result for [7 = 6* 
and V — 3.155*. The dotted line is the corresponding field- 
theoretical result (f3 2 = 0.327, E M = 0.643*). The DDMRG 
spectrum for U = 4.15*, V — 2* and the field-theoretical spec- 
trum for 1 — 1/2, Em = 0.645 are given by the solid line. 



for U — 6* and V = 3.115*, although optical gap, Mott 
gap, and total spectral weight are identical for both spec- 
tra. The DDMRG result shows that there is substan- 
tial optical weight both at E opt = 0.445* and above the 
Mott, gap Em = 0.643* while, according to field the- 
oryjlj for a ratio E opt /EM ~ 0.7 the optical conductivity 
should be dominated by an excitonic peak at u> = E opt 
with very little weight in the holon-antiholon band above 
the Mott gap Em and in the exciton-exciton continuum 
above lj » 1.7E opt ~ 1.1*. 

This disagreement between field theory and DDMRG 
results is not due to a failure of either method in the ex- 
citonic regime. It has been shown that field theory and 
DDMRG calculations for excitons agree very well in the 
extended_Hubbard model with next-nearest-neighbor re- 
pulsion.lL3 The problem is that the field-theory approach 
assumes that the low-energy optical excitations are made 
of two elementary charge excitations. In the extended 
Hubbard model ([!]), however, the conditions Em *C 4* 
and V > 2t are satisfied only close to the phase bound- 
ary V c = U/2, where low-energy excitations are CDW 
droplets made of many elementary charge excitations. 
For instance, the ionicity of the 1B~ state becomes sig- 
nificantly larger than 1 as V increases above 2* as seen in 
Table |. Therefore, the field-theory approach is not ap- 
plicable to the lattice model ([!]) with V > 2t even in the 
limit of a small Mott gap. Moreover, the extended Hub- 
bard model (Q) cannot describe a Mott insulator with a 
small gap and an exciton in the optical spectrum for any 
parameters U and V. 



V. CONCLUSION 

I have investigated the linear (one-photon) optical ex- 
citations of a one-dimensional Mott insulator, the half- 
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FIG. 12: Schematic representation of the different regions in 
the (U, V) parameter space where a particular optical excita- 
tion dominates the low-energy optical spectrum of the Mott 
insulating phase. The solid line is the boundary between the 
CDW and Mott phases. 



filled extended Hubbard model, using DMRG methods. 
Four types of optically excited states have been found: 
pairs of free (unbound) charge excitation, excitons, exci- 
tonic strings, and CDW droplets. Correspondingly, there 
are four different regimes in the model parameter space 
(U, V) depending on the nature of the low-energy opti- 
cal excitations. They are shown in the schematic "phase 



diagram" of Fig. |12[ Note that only the V = 2t line 
separating the regime of free excitations from that of 
bound excitations represents a sharp transition. The 
other dashed lines represent smooth crossover from one 
regime to another. In each regime one observes optical 
spectra with distinct features. In all cases, optical exci- 
tations are made of an even number of elementary exci- 
tations carrying opposite charges in the lower and upper 
Hubbard bands. The different types of excitations and 
optical spectra found in this system result from the resid- 
ual interactions between these elementary excitations. 

For V < 2t the low-energy optical excitations are made 
of two unbound elementary charge excitations. They give 
rise to a single continuous band in the optical spectrum, 
starting at the Mott gap Em, which thus equals the opti- 
cal gap E op t ■ The optical conductivity vanishes smoothly 
as y/iv — E opt at the threshold E opt , except for V = 2t, 
where it diverges as - E opt . For V > 2t (=> U > 

At) but V < U/3 + 0(t) and U - 2V > t, the lowest opti- 
cal excitation is an exciton (a neutral excitation made of 
two bound elementary charge excitations) at an energy 
-Eopt lower than the Mott gap Em- This exciton gives 
rise to an isolated <5-peak in the optical conductivity at 
E opt ■ Furthermore, one still finds a continuous band 



to 



starting at the Mott gap Em due to free charge excita- 
tions. For V > U/3 + 0(t) and V > 2t but U - 2V > t 
(=> U > 12t) the low-energy optical excitations are exci- 
tonic strings consisting of n cxc > 1 excitons with energy 
E(n exc ) ~ E opt + (n cxc — !)({/— 2V). The spectrum con- 



sists of several isolated 5-peaks at ui = E(n cxc ) < Em- 
In this regime one still observes a very weak band of free 
charge excitations starting at Em in the spectrum. Fi- 
nally, close to the CDW phase boundary (U — 2V < t) 
for V > 2t the low-energy excitations are CDW droplets 
and give rise to a broad band starting below the Mott 
gap E M - 

As long as optical excitations are made of a pair of 
(bound or unbound) elementary charge excitations (i.e, 
the excitation ionicity is / < 1), the optical spectra calcu- 
lated numerically with DDMRG agree perfectly with the 
analytical results obtained with a strong-coupling anal- 
ysis or with field-theoretical methods. This agreement 
confirms the accuracy and the power of the DDMRG 
method for calculating dynamical spectra in the thermo- 
dynamic limit. It also confirms the wide range of validity 
of both analytical approaches. 

Some results presented here suggest further investiga- 
tions. First, in the extended Hubbard model (|l|) no ex- 
citon exists in the regime of a small Mott gap (Em t) , 
which is relevantJor some real materials such as conju- 
gated polymers.Era It is believed-4,hatAn electron-electron 
interaction, with a longer rangal3'ta't3 or a lattice dimer- 
izatior£3o can lead to the formation of excitons in sys- 
tems with small gaps. However, the precise nature of the. 
optical excitations in such systems is still controversially 
The approach used here for the extended Hubbard model 
will enable us to determine the optical properties of these 
systems reliably. Second, excitonic strings appear in the 
linear optical conductivity spectrum for strong interac- 
tion U,V ^S> t because of the weak hybridization of ex- 
citonic strings with different S1ZCS Tioxc ■ 

Experimentally, 



excitonic strings have 
optical absorption only 



sen observed in the non-linear 
] It would be desirable to check 
if excitonic strings can be found in the linear optical ab- 
sorption of materials which are believed to be large-gap 
one-dimcnsip|Oa|l Mott insulators such as Cu oxides and 
Ni halidcsBEHlZl Last, there is a clear boundary between 
free and bound excitations in the low-energy optical spec- 
trum at V = 2t. As discussed in Ref. ^3] the nature of the 
low-energy charge excitation seems to be correlated with 
the order of the transition from the Mott insulating phase 
to the CDW insulating phase. It is likely that the tricrit- 
ical point where the transition changes from continuous 
to first order is located precisely on the line V — 2t. This 
suggests the existence of a hidden symmetry in the charge 
sector of the extended Hubbard model ([!]) at V — 2t. It 
would be interesting to investigate this feature further. 
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APPENDIX A 

In the original implementation of the charge-conjuga=j 
tion and spin- flip symmetries for DMRG calculations ,E3 
an explicit matrix representation of the superblock Ham- 
iltonian is built. This matrix can be projected onto 
a symmetry subspace with chosen parities P c and P s , 
which allows one to compute eigenstates of this symme- 
try and reduces the computer memory and CPU time 
required. In an efficient implementation of DMRG, how- 
ever, an explicit representation of the superblock Hamil- 
tonian should not be constructed (at least for quasi-one- 
dimensional systems with only short-range interactions). 
A representation in terms of tensor products aLmatrices 
uses much less memory and is also much faster .Ej Project- 
ing this representation onto a symmetry subspace slows 
down the program considerably. Therefore, insteacLof a 
projection, I use an exact diagonalization techniquet3 to 
shift the chosen symmetry subspace to lower energy. 

Let P c and P s be the charge-conjugation and spin- flip 



operators for the full lattice with eigenvalues P c = ±1 
and P s = ± 1 . As P c and P s commute with the Hamilto- 
nian H, the operator 

H' = H- X C P C - \ S P S (Al) 

has the same eigenstates as H but its eigenvalues are 
shifted, E' n — E n ± A c ± A s , where the signs ± are given 
by the eigenstate parities P c and P s . It is obvious that 
the lowest eigenstates of H' lies in the symmetry subspace 
withP c = A c / 1 A c | andP s = A S /|A S | provided |A C | and |A S | 
are large enough. Therefore, one can simply apply the 
usual ground-state DMRG approach to the Hamiltonian 
H' with appropriate values of A c and X s to obtain the 
lowest eigenstates in any symmetry sector. A similar 
approach has already been used to shift states with high 
total spin S to higher energy in a DMRG calculation.^ 
Using the method proposed recently for including a non- 
abelian symmetry group ip_.a DMRG calculation would 
be a further improvement .L3 
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